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Abstract 

An analytical expression for square- well fluid direct correlation function (DCF) obtained recently 
by Tang (Y.Tang, J. Chem. Phys. 127, 164504 (2007)) in the first-order mean spherical approx- 
imation is extended for wider well widths (2 < A < 3). Theoretically obtained direct correlation 
functions and radial distribution functions for square- well fluid with A = 2.1 and A = 2.5 are 
compared with corresponding results of Monte-Carlo simulation. 
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I. INTRODUCTION 



Recently, the systems of particles interacting with discrete potentials gained much atten- 
tion from the scientific community. Such an increased interest to this class of model systems 
is associated primarely with their ability to mimic the bulk properties of variety of complex 
fluid systems, like associating fluids, colloids, cluster particles etc. The discrete hard-core 
potentials composed of an attractive square well and a repulsive square barrier were shown 
to induce the liquid-liquid phase transition in a single- component model fluid^^i. Simi- 
lar phase transitions take place in the real substances like water, carbon and phosphorus. 
Moreover, the effective potentials with square-well plus square-barrier components may be 
used to mimic colloidal particles interactions, which according to the DLVO theory are char- 
acterized by a long-range repulsive barrier and one or possibly two attractive wells. Finally, 
inhomogeneous fluids with descrete repulsive and attractive interaction potentials represent 
another intriquing area of research because they often serve as a benchmark to study a 
variety of interesting problems such as interfacial phenomena, surface adsorption, wetting, 
capillary condensation, etc. 

The most suitable method to handle these systems seems to be the classical density- 
functional theoriy (DFT). There are several versions of the DFT derived usualy from either 
weighted density approximation or perturbation expansions. Recently, Tang and Wu^ 1 ^ 
have proposed a new DFT in which fundamental measure theory is combined with the so- 
called first-order mean-spherical approximation (FMSA). Its implementation requires knowl- 
edge of the thermodynamics and the direct correlation function of the model bulk fluid. 
In comparison with alternative approaches, DFT/FMSA possesses a number of advantages 
due to the both accuracy and simplicity of the FMSA solution. In combination with the 
DFT, the FMSA was already applied to study Lennard-Jones^, Yukawa 7 , Sutherlad^ fluid 
systems. Regarding the systems with descrete potentials, quite recently 1 - the DFT/FMSA 
approach has been applied to study the structural properties of the square-well (SW) fluid 
model defined by the pair interaction potential, 
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r < 1, 



r > A. 



1 < r < A, 
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with a unit hard-core diameter, the attractive strength e and the well width parameter 

1 < A < 2; the FMSA solution for the direct correlation function (DCF) of this model 
system has been obtained in a separate studyiS. It has been shown that DFT/FMSA gives 
good performance for entire considered range of attraction parameter 1 < A < 2. 

The fact that FMSA is linear with respect to the interaction potential allows one to use 
linear combinations of the different potential functions for which FMSA solution is known. 
Regarding the case of the SW model this means that one can employ the linear combinations 
of the square wells (s > 0) and/or square barriers (e < 0) in order to form all variety of the 
discrete potential models. Unfortunately, majority of potential applications of the discrete 
SW-based potential models (e.g., see Refs.-^ 1 ^) have the range of interaction that exceeds 
one particle diameter, that requires to consider A > 2 in Eq.(l). This means that existing 
FMSA solution obtained by Tang^ is not sufficient to have FMSA be applied in the studies 
of the complex discrete potential models. The aim of present work is to extend the recent 
FMSA solution^ for the SW model fluid with 1 < A < 2 to the case of the SW fluids with 

2 < A < 3. Obtaining of a such FMSA solution for the attractive SW model is crucial for 
the future application of both the FMSA as well as the DFT/FMSA to the systems with 
a combined (attractive SW plus repulsive SB) discrete potential. This will be shown in a 
forthcoming paper []. The actual paper is organized as follows. In the next Section II we 
outline the necessary details of the FMSA solution for SW potential with 1 < A < 3. The 
FMSA results for the direct and pair correlation function will be presented and compared 
against computer simulation data in Section III. We conclude the paper in Section IV. 

II. FMSA FOR SQUARE- WELL MODEL FLUIDS 

The general FMSA formalism, developed by Tang and Lu, is presented in their paper—. 
These authors solved the Ornstein-Zernike (OZ) equation, 

h(k) =c{k)+ph{k)c{k) 

for a one-component system of particles/moleculs of a number density p by employing the 
perturbative expansion for the total and direct correlation functions, h and c, respectively, 
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over the pair interaction energy parameter e, 

h(k) = ho (k) + eh\ (k) + ... (2) 

c(k) = c (k) + eci (k) + . . . , (3) 

where the subscript denotes the contribution of a hard-sphere (HS) reference system while 

the subscript 1 stands for the first-order perturbation term. Here and in what follows, 

in accordance with the previous notation^ all symbols with the tilde denote the three- 
dimensional Fourier transforms, 



An 
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h(k) — — I sin (kr) rh (r) dr, (4) 
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c W = ~r / s ^ n (^ r ) rc ( r ) dr, (5) 
^ Jo 

while all symbols with the hats denote the one-dimensional Fourier transforms or the Laplace 
transforms. 

By employing the Hilbert transform, the general solution of the first order OZ equation 
can be obtained. The Fourier transform of the first order contribution to the total correlation 
function h reads, 

% x {k) = ^-, (6) 



where 



P (ik) = Ul{k) - e —T Uliy)eiy dy (7) 
2Ql(-ik) 2m 7-00 (y - fc) ^ (-iy) ^ 

and Qo (ik) is the Baxter hard-sphere factorization function with a Laplace transform given 

by 

h S (s) + mjL (s) e~ s 

Qo (s) = -2— , (8) 

(1 — rj) s 6 

with rj = \kp, and 

S (t) = (1 - nf t 3 -h 67/(1 — n) t 2 + 18n 2 t - Yin (1 + 2n) , (9) 
L{t) = (l + |) t + 1 + 277. (10) 

The function U\ (k) in Eq. (0) is defined as 



00 

ikr 



U^k) = / rci (r) e~ lkr dr, (11) 
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where, according to the FMSA closure, 

c i ( r ) = —fl u ( r ) ? for r > 1, (12) 

with f3 = 1/ksT and u (r) being the pair potential. 

For smooth pair potentials u (r), for which U\ {k) ~ e~ lk (e.g., for the Yukava potential), 
the integration contour in the right-hand side of Eq. ([7]) can be closed in the upper complex 
half-plane and evaluation of the integral does not require to calculate the residues at zeroes 
of the function Q (—ik). This important simplification does not hold for the pair potentials 
that vanish for r > A (A > 1), e.g. for the SW potential. For such a potential Tang and Lu 
expanded 1/Qo i~^) 2 in a following way, 



1 {l-^fi-ikf 



Q (-ik) 2 SP(-ik) 



1-2 12T7 ; ./. e tK +3 12t7 .' e 



Jk 
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s(-ik) J v s(-ik) 

J (13) 

and performed calculations according to the scheme described in the Appendix of their 
work 16 . Since only the first term of the expansion (jTBl survives for A < 2, therefore 

Ql(ik) \(y-k)S? (-!»)/' 

where 

poo 

W x (k) = J r Cl (r)e- ikr dr (15) 

should be calculated for c\ (r) smoothly extended from r e [1, A) to r G [A, oo). For the SW 
potential we gelr^ 

W X (k) = (5e l -±^e^\ (16) 

{IK) 

Expression ( fT4l) is valid only for the pair potentials u (r) with the attractive well width 
that does not exceed one particle diameter, i.e., for A < 2. In order to extend the actual 
FMSA scheme for larger values of parameter A, the contribution due to the second term in 
the expansion ([TBI should be taken into account. After doing that we obtain the following 
result, 



2 < A < 3, (17) 



Qliik) {(y-tySH-iy) 

which may be viewed as Eq. (fT4l) supplemented with an additional term, that contributes 
to r > 2 region in the r-space. 



The thermodynamic and structural properties of the square-well fluid for A < 2 within 
the FMSA approximation have been reported in Refs.— liMS. Here, we would like to pay 
attention to the caase with A > 2, i.e., to the development of an additional term that already 
appears in Eq. ffTTl) . Let us denote it as 

h? (fc) = m ~ e "\ e s { M (-*) >l g 1 

Q5(ifc) \ (»-*)S 3 (-iy) J 

After calculating the sum of residues in above expression, we obtain the following Laplace 
transform 
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OS 00 I S»(-s) 

G 2 (U) , d^) G (ti) 



where C (s) = As 5 — s 4 , while tj correspond to the roots of S (t) = 0, and functions Gi (s) 
are, 

^.1 = -^, (20) 



4C(s)L(s)] 



ds 



+ (A - 2) C (s) L ( S ) j c {s) L {s) ^ {s) eS(A _ 2) 



Gi(s) = 2 " sm w — ' (21) 

f WWl + 2 (A - 2) d[c{ f + (A - 2) 2 C (s) L (s) ) e s ( A - 2 ) 
G 2 (.) - -- f ^ 1 (22) 

*2<gM + (A - 2) C {s) L (s) \S 2 (s) c(s)L(s)S,(s) ^ 



„ k ^2+1^-^ - - C (s) L (s) £3 (s) e s 



-3 



(a) Sf (s) 

C (s) L (s) SI (s) e< x -V 



SI (s) 

with S n (s) = d n S (s) /ds n . 

Then the Laplace transform of the first order contribution to the total correlation function 
h in the case of 2 < A < 3 can be written as, 

In the above we introduced two functions 

C (t) e^- 1 )* 



p« ( s ) = elites 



(a + t) S 2 (t) 



(23) 
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and 



(s) = e - s Res 



C (t) L (s) e( A - 2 )' 



(24) 



(s + t) S 3 (t) 

The first one was examined in more details in Ref.— , whereas the second function, in agree- 
ment with Eq. ( fT9l . reads 

C(-s)L(-s)e- s ^ 



2 ^ V s + t 

i=0 x 



53 (-*) 

G 2 (ti) G 1 (U) Go (*0 



S + t 




(25) 



Following by Tang and Lu^, the Laplace transform of the DCF can be obtained from 
the relation 



[0,oo] 



Qp(s) 



[0,oo] 



(26) 



c (s) = (3e (1 - rj) 11 1 [Q p (s) 

where superscript [0, oo] denotes the part that is nonzero only in [0, oo] region of the r space, 
and 

4 (s) = Ql (s) [p« ( s ) + 24^ 2 ) (s)] . (27) 

In order to invert c(s) to the r— space, we should ignore all [0, oo] superscripts in relation 
( [2S1 and invert it to r— space. Since we are interested only in r > region, we assume that 
c(r) = for r < 0. Note, that the contribution from the first term p^ (s) in Eq. ( 1271) was 
already considered in Ref.— . Let us examine now the second term of Eq. ([2] 



(2) 



Q 9 W (s)=24r ] QU-s)p^ (s) . 



(28) 



Taking into account that Qq (s) is analytical function in the whole complex plane, we may 
simplify the expression for p^ (s) by subtracting some analytical function. Indeed, the 
subtracted analytical function being multiplied by (— s) will result in a new analytical 
function, that has no impact on the resulting expressions for c\ (r) in the r— space. In such 
a way we replace the first term —C (— s) L (— s) e~ s( - A_2 )/S' 3 (— s) in Eq. ( |25l) by expression 

e -s(A-2)^ w h ere coefficients C\ (U), C 2 (U) and C3 (U) are 



( C 3 (tj) , C 2 (t,) , Ci(ti) 



found from the equality of residues at all poles t; of the expression 



Res 

{U} 



C 3 (U) , C 2 (^) 



s + ^ 



+ 



[S + ti 



{s + U 





= Res 




{u} 



C(-s)L\ 



■s) e 



-s{\-2) 



(29) 
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The new simplified function p'( 2 ' (s) reads 



p'lV (,,) 



E 

8=0 



C 3 (U) , C 2 (ti) d(^) 



+ 



+ 



-s(A-l) 



i=0 



(30) 



where 



d (to 

d 

^3 (ti) 



C(t,)L(t,) 
Sf(U) ' 

5? (U) + 2 ft) 



&[C(s)L(s)} s=u 3 j;[C(s)L(s)] B=ti S 2 (t i ) 
2Sf (U) 2 Sf (U) 

C (t t ) L (ti) S 3 (U) 3C (U) L (U) S 2 2 (U) 
2 



2Si(ti) 2 SUU) 

It is easy to find that the functions d (s) from Eqs. (21)-(23) can be expressed in terms 
of Cj (s) in the following way 

Go (a) = -2d (a) e^ 5 , 

d (a) = -2 [d (a) - (A - 2) d (s)] e^, 

d (a) = -2 



d (a) - (A - 2) d (a) + ^-^-d (s) 



,(A-2> 



Then the inverse Laplace transform of Eq. (1301 reads 



p'(2) ( r) 



1 2 / 

2EI G a (t,) + (r - 1) d 

i=0 ^ 



+ (r Co (^) ] e^ (r - 1} [F (r - 1) - H (r - A + 1)] 



(31) 



which is nonzero inside [1, A — 1], and H (r) is the Heaviside step function. Substituting 
(13T?j) and © into J2BD, we obtain 



n (2) n 

W _ £'(2) / n 



+ 



24r? 

d (ti) d (t 4 ) 



6 6 6 

eo,i , e i,« _ gS + e 2,» c2s 



i=l 



i=2 



i=4 



E 

8=0 



C 3 (U) 

S + U 



+ 



: s + t 4 ) 2 (s + t) 3 



e -(A-D- + 



8=1 



8=2 



+ E 

8=4 



\ 4 



E 

8=0 



d(ti) , d(ti) G (U) 



s + u {s + ur {s + uy 



(32) 
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and 



Q 



(2) 



2Ar] 



y^ + y ^ e - s + y ^ e - 2s 

i=l i=2 2=4 



E 

i=0 



C 3 (U) 
s - U 



C 2 (U) C X {U) 



(s-t 

+ y ^ 



i=4 



s — t 

2 

E 

i=0 



i=l 



i=2 



eu e - s 



G 2 {U) G X {U) G (U) 



s — 



u (s - u) (s - uy 



(33) 



where the coefficients e n j of the expansion of the Baxter function Q (s) are given in the 
Appendix. 

In order to obtain the DCF in r— space we will use the method that is similar to that 
outlined in Ref.— . Namely, from Eqs. (|32|) and (|33|) it is evident, that the key expressions 
to be inverted are the functions of following form 



s n (s + z) 1 



E -fc- + E j£=, % ± o, 

k=l k=l 

1 z = 0. 



(34) 
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with coefficients 



_1 " 



(m + n — k — 1)! 
(n - fc)! (m - i)\ z m+n-k- 
(m + n — k — 1)\ 



(m-k)\(n- \)\z m+n - k ' 
The inverse Laplace transform of these functions reads 



(35) 
(36) 



w (r, z, n, m) 



E a k 



+ E&; 



fan r k 1 g— 2r 



'* (fe-1)! ^ ^ fe (fe-1)! 
,fc=l fe=l 

\H(r), 



H(r), z^O, 
z = 0, 



(37) 



(m+n— 1)!~ 

where if (r) is the Heaviside step function. With this in hands, we may now write down in 
terms of w (r, z, n, m) the expressions for the inverse Laplace transforms of Eqs. (I32j) and 
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33]) for r > 0, 

Q? ] (r) 



6 2 



eo.i 



C 3 (tj) ^ (r - A + l,^,*,!) 



+G 2 fo) w (r - A + 1, t 3 -, z, 2) + d (t,) w (r - A + 1, t i; i, 3) H (r - A + 1) 



+££ 



G 3 (tj) w(r-X + 2, tj,i, 1) + G 2 (tj) w (r - A + 2, t,-, i, 2) 

; — i i 

i=2 j=0 

+Ci (tj) w (r - A + 2, tj, 2,3) Jf (r — A + 2) 

6 2 

&2A 



-!)• 



G 3 (tj) w (r - A + 3, t,-, z, 1) 

I — I I L 

i=4 j=0 

+G 2 (tj) w (r - A + 3, t j} i, 2) + d {tj) w (r - A + 3, t j; i, 3) if (r — A + 3) 

6 2 

c ,-, ■ r 

G 2 (tj)w(r - l,tj,i,l) 



eo,-, 



i=i j=o v 



+d (i,) w(r-l, tj, i, 2) + G (tj) w{r-l, tj, i, 3) H (r - 1) 



6 2 



e i,-. 



G 2 {tj)w(r, tj,i, 1) 



8=2 j=0 ( 1) 

+d fo) iw (r, 2) + G (tj) w (r, t„ i, 3) if (r) 



1 6 2 



G 2 (tj)w(r + l,tj,i,l) 



/-:] j=0 ( ' 

+Gi (tj) w (r + 1, tj, i, 2) + G {tj) w(r + l, tj, i, 3) H (r + 1) , 



(38) 
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gg (-r) 
24ry 



6 2 



p(2) (_ r ) _ ^ ^ e0i . c 3 u»(r + A- 1, z, 1) 



i=l i=0 

-G 2 fa) w (r + A - 1, -tj, i, 2) + d to (r + A — 1, -t 3 -, i, 3) if (r + A - 1) 

6 2 

- X! X) ei >* [ ^ w ( r + A - 2, i, 1) 

i=2 j=0 

-C 2 (tj) w(r + \-2, -tj, i, 2) + d (tj) w(r + \-2, -tj, i, 3) if (r + A - 2) 

6 2 

-XX e2 >* ^3 ^ (»" + A -3, -tj,i,l) 

i=4 J=0 

-C 2 (t 3 -) iu (r + A - 3, -t i9 1, 2) + d (tj) w (r + A - 3, -t,-, z, 3) H (r + A - 3) 

1 6 2 r 
_ 2^5Z e °'' G 2(^)w(r + 1, -tj,i, 1) 

i=l j=0 

-G x (tj) w(r + l, -tj, i, 2) + G (tj) w(r+l, -tj, i, 3) H (r + 1) 

6 2 



- g X X ^ [ G2 w ( r ' *' X ) 

i=2 j=0 

—Gi (tj) w (r, -tj,i, 2) + G (tj) w (r, —tj,i, 3) H (r) 

6 2 



- g X Yl e ^ [ G2 w ( r ~ 1; *' ^ 

i=4 jf=0 

-Gi (tj) w (r - 1, i, 2) + G (tj)w(r- l,-tj,i,3) H (r - 1) . (39) 



The above expressions are then substituted into inverse of Eq. (127)) . 



Q P (r) = Q?> (r) + Qf (r) 



(40) 



where Qp 1 ^ (r) is the inverese Laplace transform of the function Q p ^ (s) = 2Ar]Ql (—s) (s) 
and describes the contribution of the first term of the expansion (ITTj) ; the expression for this 
contribution was given in Ref.— . Then the first-order perturbation contribution to the DCF 
of the SW model fluid with the radius of interaction that exceeds two particle diameters 
reads, 

rc x (r)=(3e(l- V ) 4 [Q p ( r )-Q p (-r)\. 

The equations (1381) and ( 1391) have been presented here in a similar manner as it has been done 
in the work of Tang and Lu^ for the case of the SW model fluid with the range of interaction 
that does not exceed two particle diameters, i.e. for A < 2. The resulting expression for 
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c\ (r) in the case of A > 2 is consequently rather long, but also quite straightforward. It 
may be summarized as follows, 



rci (r) 



0. 

per, 



r > A, 
1 < r < A, 



(41) 



f3e{l-r ) f[Q p {r)-Q p {-T)\, r < 1, 



where Q p (r) is given by Eq. (j4T)l) . The direct correlation function is discontinuous at r = 1 
and r = A. The continuity of the DCF inside the core, r £ [0, 1], is evident from Eqs. fl32|) 
and (13"3"j) . where the terms, that might contribute to the discontinuities at r = A — 2 and 
r = 3 — A, are at least of the second order in 1/s. 



III. RESULTS AND DISCUSSIONS 

To examine the evolution experienced by the DCF of the SW fluid upon increase of the 
attractive well width, in Fig.[TJwe show the first-order DCF c\ (r) evaluated for nine different 
values of width parameter A in the range from A = 1.5 to A = 3.0. All calculations are 
performed for fixed density p = 0.75 and temperature j3e = 0.5. Note, that the dependence 
of Ci (r) on the attractive well width A in the core region r < 1 is not trivial. 

It is very well illustrated by the behavior of the value of c\ (r = 0) . One can see that 
value of Ci (r = 0) is first decreasing upon A increases from A = 1.5 to A ~ 1.7 and after 
this starts to grow with A increases. However, after reaching the value of A ~ 2.2, c\ (r = 0) 
starts to drop till the value of A ~ 2.6 after which it begins to grow again. 

In Figures [2] and [3] we show the total DCF c (r) and total radial distribution function 
(RDF) g (r) = h x (r) + 1 of the SW fluid with attractive well width A = 2.1 and 2.5 at 
the density p = 0.75. Both functions are obtained within the FMSA and from Monte- 
Carlo simulations. The FMSA results for the first-order correction term hi (r) to the RDF 
were obtained numerically. The correspondin DCF and RDF of the hard-sphere reference 
system, cq (r) and go (r), respectively, were obtained from the first-order GMSA theory of 
Tang and Lu^. The results for the intermediate and low densities are not shown, since at 
such conditions the considered SW fluid seems to be unstable. Good agreement between 
Monte Carlo simulation data and FMSA theory results could be observed for the SW fluid 
with A = 2.5 (see Fig. [2]). The same doesn't hold for the SW fluid with A = 2.1. The 

12 



discrepancies between FMSA and simulation predictions are clearly visible for both DCF 
and RDF at the particle contact r = 1 and at the point of discontinuity at r = A = 2.1. 
These shortcomings in the case of RDF can be corrected by employing exponential (EXP) 23 

9 EXP (r)=g (r)e h ^ (42) 

or linearized exponential (LEXP)^ 

g LEXP (r)=g (r)(l + h 1 (r)) (43) 

approximations. Indeed, the LEXP approximation significantly improves the RDF of the 
SW fluid, giving better contact value at r = 1 and improving RDF values on both sides of 
the discontinuity at r = A = 2.1 (see Fig. [3], right panel). 

In Fig. H] we present the inverse reduced isothermal compressibilities of the SW fluids 
that was obtained from 

— = (3 ( ^ ] = 1 - Anp [ c (r) r 2 dr, (44) 
Xt \°P J T J 

The DCFs of the reference system of hard spheres were obtained from the first-order 
GMSA theory of Tang and Lu 2 ^. From Fig. H] we note that at the high densities (rj « 0.5) 
the compressibility of the SW fluid with A = 2.5 increases rapidly and almost reaches the 
level of the fluid A = 2.0. Mathematically, this can be explained by positive contribution into 
compressibility of the first-order DCF in the r < 1 region, which compensates the negative 
contribution due to the potential well. 

On Figj5] we present spinodal curves and critical points for few long-range SW fluids 
obtained from condition — and (jHJ). 

IV. CONCLUSIONS 

We extended the FMSA theory to deal with the square-well fluids of large wells width 
(1 < A < 3) than those studied previously and obtained analytical expressions for the DCF. 
We obtained reasonable agreement of the DCF and RDFs with MC simulation data and 
proposed exponential (EXP) and linearized exponential (LEXP) approximations to correct 
RDF values close to the contact. 
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APPENDIX A: EXPANSIONS OF THE LAPLACE TRANSFORM OF THE 
BAXTER FACTORIZATION FUNCTION 

Laplace transforms of the Baxter hard-sphere factorization function and square of the 
transform are 

3 3 

4)(*) = i + E- + E- e ~ s > ( A1 ) 



and 



where 



and 



i=i Si i=2 Si 

i=l i=2 i=A 



_ 12r, (1 + jj) _ 12,) (1 + 2,)) 
h 

a 1 = b 2 - —, a 2 = b 3 - b 2 , a 3 = -b 3 , (A3) 



e ,i = 2ai, e ,2 = a\ + 2a 2 , e ,3 = 2a 3 + 1a x a 2 , 

e , 4 = 1a x a 3 + a 2 ., e , 5 = 2a 2 a 3 , e 0)6 = Og, 

ei,2 = 26 2 , ei j3 = 26 3 + 2ai& 2 , e M = 2ai& 3 + 2a 2 6 2 , 

ei,5 = 2a 2 6 3 + 2a 3 6 2 , ei )6 = 2a 3 6 3 , 

e 2 , 4 = &2> e 2,5 = 26 2 6 3 , e 2i6 = 63. 
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FIGURE CAPTIONS 



Fig. [T] First-order DCF for several A values with e/3 = 0.5 and p = 0.75. 

Fig. [1 Full DCF (left panel) and full RDF (right panel) for SW fluid with 
A = 2.5, e(3 = 0.5 and p = 0.75. Lines denote FMSA results, symbols denote 
simulation results. 

Fig. M Full DCF (left panel) and full RDF (right panel) for SW fluid with 
A = 2.1, e(3 = 0.5 and p = 0.75. Lines denote FMSA results, symbols denote 
simulation results. Dotted line denotes linearized exponential approximation 
(LEXP). 



Fig. 3] Inverse reduced isothermal compressibility x r x = (3 of the SW 

fluids at e/3 = 0.5 and e(3 = 0.25 . 

Fig. m Spinodal curves for SW fluids with A = 2.0, 2.2, 2.4, 2.5, 2.6, 2.8, 3.0. 
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